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Abstract. We compute the four loop term of the mass anomalous dimension in the two dimen- 
sional Gross-Neveu model in the MS scheme. The absence of multiplicative renormalizability 
which results when using dimensional regularization means that the effect of the evanescent 
operator, which first appears at three loops in the 4-point Green's function, has to be properly 
treated in the construction of the renormalization group function. We repeat the calculation of 
the three loop MS /3-function and construct the /9-function of the evanescent operator coupling 
which corrects earlier computations. 
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1 Introduction. 



The Gross-Neveu model is a two dimensional asymptotically free renormalizable quantum field 
theory whose basic interaction is a simple quartic fermion self-interaction, It has been 
widely studied since its introduction in [T], as it has many interesting properties which are 
readily accessible given the space-time dimension the model is defined in. For instance, unlike 
the same interaction in four dimensions it is renormalizable in two dimensions and the generation 
of mass dynamically has been observed and studied in the large expansion, [1]. Moreover, one 
property of interest is that it possesses an S'-matrix whose exact form is known, [21 [3], whence 
the mass gap is known exactly, [4j, in terms of the basic mass scale of the theory, A^^. Aside 
from these features the model itself underpins several problems in condensed matter physics. 
For instance, in the replica limit it is equivalent at the critical point to the two dimensional 
random bond Ising model. (See, for example, the review [5].) Necessary to study the fixed point 
properties for such physical problems is knowledge of the renormalization group functions in 
some renormalization scheme, such as MS. These have been computed to three loops in MS 
over a period of years. In [1] the one loop /?- function was computed demonstrating asymptotic 
freedom. This was extended to two loops in [6], whilst the three loop /3-function appeared more 
or less simultaneously in ^ and Though the method of computation in both articles was 
significantly different. For instance, given the quartic nature of the sole interaction it can be 
rewritten in terms of an auxiliary field producing a trivalent interaction with the introduction 
of an auxiliary field. This was the version of the theory used in [8j, as well as at two loops in 
[6], not only to deduce the /3-function but also to study the effective potential of the auxiliary 
field at three loops. However, renormalization effects will generate a quartic interaction. So 
[8] had in effect to handle the intricate problem of renormalizing a version of the theory with 
two independent couplings. The three loop MS /^-function of the original theory was eventually 
extracted when the effect of the newly generated interaction was properly accounted for in the 
renormalization group equations. By contrast, in [7j the purely quartic version of the theory was 
treated with a massive fermion. The agreement of both three loop results was a reassuring non- 
trivial check on the final expression. At four loops only the wave function renormalization has 
been computed in [9j and later verified in [lOj. Although apparently one loop further than the 
/3-function of [7l[8], or mass anomalous dimension, [TT], the fermionic nature of the interaction 
means that the anomalous dimension begins at two loops since the one loop snail graph is zero 
in the wave function channel of the 2-point function. Thus in effect the four loop wave function 
is a computation on a par with the three loop mass anomalous dimension and /3-function. 

Given the range of problems which the Gross-Neveu model underlies, it is the purpose of 
this article to start the programme of completing our knowledge of the four loop structure by 
computing the mass anomalous dimension at this order in the MS scheme. This may appear to 
be the least important of the two outstanding quantities. However, as will become apparent from 
the calculation we will describe, the nature of the model is such that there are several difficult 
technical issues to be dealt with en route which do not arise in other four loop calculations 
in other important theories. Therefore, in one sense we are testing the viability of computing 
renormalization constants at four loops in an example in the Gross-Neveu model which contains 
a moderate number of Feynman diagrams rather than the 1000 plus graphs which will occur 
in the 4-point function renormalization and with a sensible investment of time. Moreover, the 
experience gained in such an exercise will be invaluable for any future such coupling constant 
renormalization. It will transpire that we will require properties of Feynman integrals in two 
space-time dimensions higher than the one we compute in, in order to determine the final 
renormalization group function. Whilst this may not be the unique way to determine this, 
it will suggest the importance of the structure of higher dimensional Feynman diagrams in 
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complementing practical lower dimensional calculations and potentially equally vice versa at 
high loop order. Our main tool of computation will be the use of dimensional regularization 
in d = 2 — e dimensions where the relevant part of the 2-point function is written in terms of 
basic massive vacuum bubble graphs. Such a calculation could only be completed with the use 
of computer algebra and invaluable in this was the symbolic manipulation language Form, |12j . 
At this level of loop order automatic Feynman diagram computation using computers is the only 
viable way of keeping a reliable account of the algebra within a reasonable time. Several other 
computer packages were also required. 

As this is the first four loop calculation which involves four terms of the renormalization 
group functions, we will have to revisit and redo the earlier calculations using the same approach 
consistently as that which will be used here at four loops. For reasons which will become evident 
later this has entailed us carrying out the full renormalization of the 4-point function in the 
theory with a massive fermion, where the mass will act as a natural infrared regulator. In [7] 
and [8] these articles centred on the derivation of the three loop /3-function itself, which unlike 
most field theories, is not the same as fully renormalizing the underlying n-point function. For 
the Gross-Neveu model this was first observed in [13^ [H] with explicit three loop calculations 
for the full 4-point function given for a massless version of the theory discussed in the series of 
articles \15 \ [T7] and examined in fL8\ for the massive version. However, neither computation 
was in complete agreement as to the final structure of the 4-point function renormalization. Prior 
to considering any four loop mass anomalous dimension this discrepancy needs to be resolved in 
one way or another which is a secondary consideration for this article. Any four loop /3-function 
computation would also require this resolution but constructing the mass anomalous dimension 
in a consistent way is an easier environment in which to check any final explanation. 

The paper is organised as follows. We review the current understanding of the renormal- 
ization properties of the Gross-Neveu model in section 2. Given this we discuss the three loop 
vacuum bubble integrals needed for the full three loop renormalization in section three. There 
we resolve the discrepancy between [17J and [18j in the 4-point renormalization in section 4. 
Section 5 extends the discussion of the relevant vacuum bubble computations to four loops with 
the four loop mass anomalous finally being derived in section 6. Concluding remarks are given 
in section 7. 

2 Preliminaries. 

We turn now to the specific properties of the Gross-Neveu model. The bare two dimensional 
Lagrangian is, [T|, 

L = ii,lc^i,l - mo%% + ^9oi4i^hf (2.1) 

where the subscript o denotes a bare quantity and g is the coupling constant which is dimen- 
sionless in two dimensions. Unlike [71 [Tl| which used the symmetry group 0{N) we take the 
SU{N) version of the theory so that the fermion field ip^ is a complex Majorana fermion with 
the former property deriving from the fields taking values in the group SU{N) with 1 < i < N . 
An advantage of the choice of the group SU{N) is that the Feynman rule of (j2.1|) involves 
two terms unlike the three of the 0{N) case. At four loops this reduces the number of terms 
needed to be substituted when the Feynman rules are implemented in Form which speeds up 
the calculation. We choose to work with the massive version where m is the mass. This is 
primarily because we will need to redo the three loop renormalization of the coupling and the 
presence of a non-zero mass will ensure that any resulting divergences are purely ultraviolet 
and not deriving from spurious infrared infinities which could arise when external momenta are 
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nullified in the basic divergent 4-point Green's function. In two dimensions the theory (12. ip is 
renormalizable to all orders in the coupling constant and is asymptotically free. Specifically we 
note that the MS scheme renormalization group functions of the model, as they currently stand 
are, P El El El El El HZ] , 

7(5) = (2iV-l)#^ - (iV-l)(2iV-l):^' 



87r2 ' 'I67r3 
+ (47V2-14Ar + 7)(2iV-l)^^ + 0(/) 

7^(5) = -(2iV-l)|- + (2iV-l)|^ + (4iV-3)(2iV-l)J^ + 0(/) 

m = {d-2)g - {N-lf- + {N-l)^ + (iV-l)(2iV-7):^' 



TT 



27r2 ^ ^167r4 
+ 0(/) (2.2) 

where 7(5), 'Jmig) and /3((7) are respectively the field and mass anomalous dimensions and 
the /3-function. Their formal definitions will be discussed later. Although several terms were 
determined for the 0{N) version of the model, we have converted the previous computations to 
the SU{N) model whence the free field case emerges when iV = ^ as indicated by the vanishing 
of ^{g) and 7m (5) for this value. 

In the computations deriving ()2.2p the main strategy was to dimensionally regularize ()2.ip 
in d-dimensions and determine the renormalization constants as poles in the deviation from 
two dimensions. Here we will take d = 2 — e with e being regarded as small. Whilst the 
correct renormalization group functions emerged at three loops, ^ overlooked a novel feature 
of the dimensionally regularized Lagrangian which was explicitly discussed in [151 [T6] after 
the observation in [13\ I14j . Basically (j2.ip ceases being multiplicatively renormalizable in d- 
dimensions but crucially retains renormalizability. This is not a property solely restricted to the 
Gross-Neveu model but is a feature of any two dimensional model with a 4-fermi interaction such 
as the abelian and non-abelian Thirring models, [18]. At a certain loop order, which is different 
for different models, evanescent operators are generated through the renormalization which 
are non-trivial in d-dimensions but which are absent or evaporate in the limit to strictly two 
dimensions which corresponds to the lifting of the regularization. A comprehensive study of this 
problem was provided for general 4-fermi theories in [131 [T4] and we recall those features which 
are relevant for our ultimate goal. The same problem in four dimensions has been considered 
in [191 120j . Though there 4-fermi operators are of course non-renormalizable and treated in the 
context of effective field theories. 

First, in d-dimensions the basis of 7-matrices based on the Clifford algebra 

{7^7'^} = 27?^^^^ (2.3) 

has to be extended to the set of objects T^^j"'^" for integer n > 0, [T3l [m [151 [HI EZ], which is 
totally antisymmetric in the Lorentz indices and defined by 

where we use the convention that the square brackets include division by n\ when all possible 
permutations of the 7-strings are written explicitly. Then T'^^^"^" form a complete closed basis 
for 7-matrices in d-dimensions where F^q) is the unit matrix. Hence one can immediately see that 
the most general multiplicatively renormalizable 4-fermi theory using dimensional regularization 
in d-dimensions is, [IS [M [13 [TSl [Hj , 

1 00 

L = - moiji^^o + ^E5(n)oV5orj;Y'''"V'^V5^r(„)^^...^>^ (2.5) 

n=0 
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where there is an infinite number of (bare) couphngs q with ^(o) = 9 identified as the original 
one of the Gross-Neveu model (|2.ip . Though the Gross-Neveu model strictly will correspond 
to the case where (7(1) = 5(2) = as Fj^-^^ and T'^^-^ are not evanescent. Given (12. 5p . there 
are several points of view depending on the problem in hand. If (j2.5p is the most general 
renormalizable theory in d-dimensions, then in principle for the Gross-Neveu model one must 
begin with ()2.5p but omit and 5(2)- This will produce renormalization group functions 
dependent, in principle, on all evanescent couplings. The true renormalization group functions 
of the original theory would eventually emerge from this multiplicatively renormalizable theory 
by setting = for n > 3 at the end, fTSj [Ml [I5l [Ml jTYj. Clearly this would involve a 
significant amount of calculation much of which would be redundant in the production of the 
final renormalization group functions. From a practical point of view there is a less laborious 
route to follow if one abandons the insistence on multiplicative renormalizability, \17\ [T8] . Then 
operators such as 

On = iV5^Ff--'''"^^^^rWMi-M>' (2-6) 
for n > 3 will be generated with (7(0) = 9 dependent coefficients. The problem for this point 
of view then becomes one of how to extract the true two dimensional renormalization group 
functions. It turns out that a formalism was developed in \13\ [T^ and used in |18[ [TO] for this 
evanescent operator issue. In essence the true renormalization group functions are not strictly 
determined from what we term the naive renormalization constants. By these we mean those 
required to render 2 and 4-point functions finite. Instead these naive renormalization group 
functions need to be amended by the effect the evanescent operators have on the divergence 
structure in d-dimensions relative to two dimensions. In [13\ IT^ such a projection formula was 
introduced which involves projection functions, p^^\g), Pm{9) and C^^\g), where the index 
k ranges over the evanescent range A; > 3. These functions quantify the effect the evanescent 
operators have on the derivation of the renormalization group functions. The derivation of the 
projection formula is given in [T3] and applied additionally in |181 110|. We recall that it is 



' 9(i)=0,rf=2 



where the normal ordering symbol, N , is included, |131 [T^ [2T| I22j . The relation strictly only 
has meaning when inserted in a 2 or 4-point Green's function. In other words one inserts the 
evanescent operator of the left side of (j2.7p in a Green's function and evaluates it using the naive 
renormalization constants to yield a finite expression. Equally one inserts the left side of (12. 7p 
into the same Green's function to the same loop order and renders it finite. Then the coefficients 
of the perturbative expansion in the coupling constant g are chosen order by order to render the 
equation consistent at each loop order after one has set d = 2. This procedure is denoted by 
the restriction {^((j) = 0, d = 2} on both sides of (j2.7p . Once the explicit projection functions 
have been determined to the appropriate order, then the true renormalization group functions 
are given by, [131 [H! ■, 

Pig) = Pig) + J2C^'\9)Pki9) 
k=3 
00 

7(5) = 7(5) + Y.P^'^(9)Pki9) 



k=3 

00 



imig) = imig) + P^^\9)Pkig) 

k=3 
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where ~ denotes the naive renormahzation group functions. 

For the Gross-Neveu model the first appearance of an evanescent operator is at three loops 
which was originally observed in |14[|17| . Whilst this postdates the three loop MS /3-functions of 
[71 E] the latter are unaffected by the generation of O3 since it occurs with a coupling dependence 
of 5^. So that coupled with C^^\g) it will only affect the /3-function itself at four loops. Equally 
the mass anomalous dimension of pAj does not feel this evanescent operator presence until four 
loops either. We refrain from quoting the value of the associated /^-function, (33(g), until later. 
This is primarily because there are two completing values given in [17J and [18j. In the former the 
renormalization was deduced in a massless version of (j2.ip where is was claimed that only ladder 
style diagrams were the origin of O3. In that instance the nullification of two external momenta 
in the associated 4-point function should not have resulted in spurious infrared singularities. 
Whilst a P3{g) was determined, it involved (^(3) which was not found in [18] which used the 
massive version, ()2.ip . where is the Riemann zeta function. This clearly avoided infrared 
singularities when all the external momenta were nullified in the 4-point function. Moreover it 
was claimed that the diagrams leading to O3 were akin to those analysed in [T^ but with no 
C(3) appearing in the published value of (33(g). Though both calculations agreed on the rational 
part of (33(g). The discrepancy between both computations needs to be resolved and a four loop 
calculation which requires (33(g) explicitly to obtain the true renormalization group function will 
provide a non-trivial forum in which to achieve this. The correct expression for (33(g) will be 
crucial for the four loop /3-function. Given this structure of the Gross-Neveu model we can now 
write down the renormalized form of (j2.ip we will use. It is, jlT^ llSj. 

L = iZ^^'^^' - mZ^Zmi^'ij' + ^gfi'ZgZl(iP'r)^ + ^gfi' Z 33 Z^ {tp'T^s^i^'f (2.9) 

where the renormalized quantities are defined from their bare counterparts by 

1 

ipo = tpZ^ , mo = mZm , go = gZg^i^ (2.10) 

in d-dimensions and Z33 absorbs the infinity associated with the generation of O3 at this order. 
Unlike Z^, Z^a and Zg its coupling constant expansion does not commence with unity. Though 
we stress that ()2.9p is valid only for 2-point calculations to four loops. Only by renormalizing 
the 4-point function at four loops would the full evanescent operator structure at that order 
emerge. For instance, it is not inconceivable given the 7-matrix structure of the four loop 4- 
point function that a O4 evanescent operator will be generated. From these renormalization 
constants the naive renormalization group functions are given by 

d ~ d 

1(9) = '"^^"'^V' ' 7m (5) = - (3(g)^\^Zm 

P(g) = (d-2)g - gP(g)-^lnZg (2.11) 

to the order we are working to. For P3(g) one deduces its explicit form from the simple pole in 
e via standard methods, [14J. Thus in the context of (j2.8p and these observations, we note that 
for the mass anomalous dimension the result of [H] for pm (g) is 

pt\9) = - ^ + O(g^) . (2.12) 
vr 

The higher terms are not required since the first term of (33(g) is O(g^). 
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3 Three loop calculations. 

We begin this section by discussing the computational strategy. To determine the mass anoma- 
lous dimension for (j2.ip we consider the 2-point function for the massive theory. In [9] the four 
loop MS anomalous dimension was calculated and independently verified in [TO]. Therefore, we 
assume that result for Z^. However, this is effectively a three orders calculation since the one 
loop snail of Figure 1 corresponding to {ipaip)^^^ {—p)) has no non-zero contributions in the ^ 
channel for the massive or massless Lagrangians where p is the external momentum. Moreover, 
since for this case one is interested only in Z^, it sufficed to consider the massless theory whence 
one only has to determine massless Feynman integrals. The component involving can be 
deduced by multiplying all diagrams by |^ and taking the spinor trace. For the mass dimen- 
sion one cannot follow this strategy. Not only because the one loop diagram contributes but 
also because its contribution to Zm requires the presence of the mass itself. Therefore unlike 
the determination of Z^ one cannot neglect the snail graph of Figure 1 at one loop as well 
as the graphs where snails appear as subgraphs at higher order. However, given that one is 
only interested in the mSj^ channel of {ipa{p)'4^^ {—p)) the Green's function can be analysed by 
nullifying the external momentum. Taking the spinor trace produces Lorentz scalar integrals 
but with tensor structure resulting from internal momenta contractions. The presence of the 
common mass m automatically protects against the appearance of spurious infrared infinities 
and relegates the problem of determining the ultraviolet structure to mapping the integrals with 
internal momenta contractions to a set of basic master vacuum bubbles at each loop order. The 
problem of studying massive vacuum bubbles in four dimensions has received much attention 
over the years, culminating in, for example, the Matad package at three loops, [23], and the 
comprehensive study by Broadhurst of all combinations of massive and massless propagators 
in the Benz or tetrahedron topology, [25]. The analogous problem relative to two dimensions 
has not been treated as systematically. Though the main results to three loops have appeared 
within various articles. Additionally, at four loops we will have to handle new integrals for 
topologies which do not simply break into products of lower loop vacuum bubbles. The main 
difficulty lies in having to handle the tensor structure emanating from the fermion propagator. 
Throughout we have made extensive use of the symbolic manipulation language Form, ^12j, in 
which to code our algorithm where the contributing Feynman diagrams to the 2 and 4-point 
functions are generated automatically with the Qgraf package, [23]. To summarize there are 
1 one loop, 2 two loop, 7 three loop and 36 four loop graphs for the 2-point function. For the 
4-point function there are 3 one loop, 18 two loop and 138 three loop graphs to renormalize. 



At three loops we now summarize the few master (massive) integrals which will be of interest 
to us, in our notation and conventions. First, Figure 2 denotes the basic vacuum bubbles to two 
loops. We define the first graph of Figure 2 as 




Figure 1: One loop contribution to the 2-point function. 




(3.1) 
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We work in Minkowski space and choose to include a factor of i with each integration measure 
which is abbreviated by 

(3.2) 



Ik J i2TTY ■ 

The integral / is trivial to deduce from the Euler /3-function as 



Hence the middle graph of Figure 2 is I^. 






Figure 2: One and two loop vacuum bubbles. 
The final graph we denote by A(0), [11], where 

A(0) = ^ f (3.4) 

and ^ 

•^^^'^ " 'Ik [k^ - m^][{k - p)^ - m^] ^^-^^ 

is the basic one loop self-energy bubble. There is a sequence of integrals related to J{p) defined 
by 



Jaf^ip) - [^2_^2]a[(^_p)2_^2]/3 M 

where we choose J2i(p) = K{p'^). In this form a Gauss relation of the hyper geometric functions 
gives the relation 

(p2_4m2)K(/) = j(o) - (d-3)J(/) (3.7) 

with 

since explicit calculations produce 

J( 2^ r(2-d/2) / W-pn ^/^-^ / d 1.3 p^ \ 



and 



, _ r(3-d/2) /W-p^y/^-' / d 1.3. ■ 

Likewise, at the subsequent level 

(p2_4m2)(j22(p2) + 2J3i(p2)) = 2K(0) - 2(d-5)E:(p2) (3.11) 
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whence 



{d-6) 
_ jd-Q) 

p2 



(p2 - 2m2) 

p2^p2 _ 4j7),2^ 



+ 



p2^p2 _ 477j2"J 



(i^(0) - (d - 5)ir(/ 
K{0) - {d-b)K{p^) 



(3.12) 



These rules are used extensively for the two and higher loop Feynman integrals. The integral 
A(0) is finite in two dimensions and can be evaluated in an expansion in powers of e as 



A(0) 



9S2 



+ 0{e) 



(3.13) 



where S2 = (2\/3/9)Cl2(27r/3) with Cl2(a;) the Clausen function. The analogous four dimensional 
vacuum bubble also contains S2 in its finite part but is divergent. In principle the 0(e) term of 
(j3.13p can be deduced. However, throughout our computations we left A(0) itself unevaluated 
since on renormalizability grounds it must be absent from the final renormalization constants at 
higher loops. This is because if the 2-point function did not have its external momenta nullified 
then the integral A(p) would emerge, where 



A(p) 



J(fc2 



(3.14) 



Clearly such a non-local function of the external momenta could not be retained when all the 
counterterms are included. 








Figure 3: Three loop vacuum bubbles. 

At three loops there are several more basic topologies. If one ignores for the moment the 
complication due to the presence of internal momenta contractions in integral numerators, then 
the basic graphs are given in Figure 3. The first involves / and its derivatives with respect to 
m?. Also the second is a variation on A(0) and we note that the two loop subgraph is 



J(A;2 



[A;2 



m 



212 



(^^-3) 
2>m? 



A(0) 



(3.15) 



which is established by differentiating A(0) with respect to m?. Aside from the Benz topology 
the remaining vacuum bubbles of Figure 3 are related to the integrals i </^(A:^)/[A;2 — m^], 
i j2(A:2) and i Jj^ J{k'^)K{k'^). Similar to A(0) these are finite in two dimensions but their 
values are required at four loops when multiplied by counterterms. As only the leading term in 
e is required in each case it transpires that we can set d = 2 and use the fact that in Euclidean 
space, denoted by the subscript E, 



JEi-k") 



d=2 



sinh^ 



J(0) 



d=2 



(3.16) 
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upon the change of variables fc^ 
instance, 



Wsinh2(6'/2) where J(0)| 



d=2 



l/(47rm2). Then, for 



27r d=2 Jq 

This can be evaluated from standard integrals to give 

7C(3) 



de 



sinh^ 



+ Oie) 



Equally we find 



647r3 



^ 11C(3) 

Though in d-dimensions one can derive the relation 

(3d - 8) 



i / J{k')K{k') 



+ 0(e) 



+ 0{e) 



J\k^ 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



The presence of K{k'^) in several of the master integrals with topology similar to those of Figure 
3 produces similar finite integrals whose finite part is required and which is determined in an 
analogous way. We note that 



J(A:2 



k [k'^ - 4m2] 

J^(fc^) 
k [A;2 - 4m2] 



ln(2) 
27r 
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J(0) 



d=2 



+ 0(6) 



C(3) - ln(2) 



J2(0) 



27r 



+ 0(6) 



(3.21) 



d=2 



For the mass anomalous dimension these basic vacuum bubbles suffice to determine the renor- 
malization constants to three loops. Given the nature of the 4-point interaction in (|2.ip the 
Benz topology does not occur in the 2-point function at this loop order. 

Having discussed the basic scalar master integrals which result we briefly note the algorithm 
dealing with the numerator structure of the integrals. This has been systematically quantified 
in [11]. However, we note that repeated use of 

1 
2 



kl 



k^ + f 



[{k-lf 



m 



m 



(3.22) 



and then 



k'^ = [k^ - w?] + , 
in each contributing topology where there are [/c^ ■ 



I' 



-m 



and [{k-l)^- 



(3.23) 

■ m^] propagators 



already. This is done in such a way that powers of kl can remain when all mixed [{k — I) — m ] 
propagators are absent and one does not then continue substituting for kl. In such integrals one 
can use Lorentz symmetry in the k and I subgraph integrals to redefine even powers of kl as 
proportional to /c^/^ or zero if there are an odd number of factors of kl. Then ()3.23p is repeated. 
Consequently several variations in the basic bubble graphs of Figure 3 emerge and we note that 

4 



i I {efj\k') 



8m^ 
3{3d - 4) 



(5d-6)I^ + 2dm^i [ J^{k'^ 
Jk 



(3.24) 



where these are exact and no finite parts have been omitted since these are crucial for the 
next loop order. In essence this summarizes the key ingredients in the algorithm for evaluating 
the three loop mass anomalous dimension which has been coded in Form and reproduces the 
previous three loop MS Gross-Neveu mass anomalous dimension. 
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4 Three loop 4-point function renormalization. 



At this point we turn to our secondary aim which is to resolve the discrepancy in the renor- 
mahzation associated with the generation of O3. This requires the complete determination of 
the 4-point function divergence structure at three loops. Whilst the algorithm to do this is 
very similar to that of the 2-point function there are several key differences. First, the 4-point 
function divergences will be independent of the external momenta which means that they can be 
immediately nullified. The mass again protects against spurious infrared divergences. However, 
we cannot now take the Lorentz traces since that would prevent one from seeing the emergence of 



anyr^3) . 



(3) /xi/i 



o- 7-matrix structure. Instead we have to retain 7-strings and also unentangle 



the internal momenta within them. Hence one decouples the Feynman diagrams into 7-strings 
and Lorentz tensor vacuum bubbles. At one and two loops the resulting tensor integrals for the 
whole integral can be straightforwardly reduced by noting that at one loop 



Ik d Jk 

where k is the sole loop momentum and at two loops 

1 



(4.1) 



kl 



d{d-l){d + 2) 

f f2{k,l) [[{d + 1)^1.^2^3-^4 - k,.k3k2.ki - ki.k^k^M r?WM2^M3M4 
Jkl 

+ [{d + l)ki.k^k2.k^ - k1.k2k3.ki - k1.kik2.k3] r,A'iA'3^/.2M4 

+ [{d + I)ki.kik2.k3 - k1.k2k3.ki - k1.k3k2.k4] ryW/^^^'^aMs] 

(4.2) 

where ki G {k, 1} and in the Lorentz tensor of the integrand all possible combinations of the two 
internal momenta are covered. The functions fii{ki}) represent the various possible propagator 
combinations. For clarity we have included the dot of the scalar products explicitly. At three 
loops the situation is complicated by the observation that the extension of both these formula 
gives 



klq 



ki^ k2^ ki^^ k^'^k^'^ kg*^ f3{k, I, q) 



(i((i- l)(d-2)(fi + 2)((i + 4) 

f3{k, l,q) {d^ + M- 2)ki.k2k3.kik^,.kQ 

klq '- 

— {d + 2)ki.k2k3.k^kQ.ki — [d + 2)ki.k2k3.kQki.k^ 

— {d + 2)ki.k3k2.kik^.kQ + 2k1.k3k2.k5kQ.ki 
+ 2k1.k3k2.kQki.k5 — {d + 2)ki.kik2.k3k5.kQ 
+ 2k1.kik2.k5kQ.k3 + 2k1.kik2.kQk3.k5 

+ 2k1.k5k2.k3ki.kQ + 2k1.k5k2.kikQ.k3 

— {d + 2)ki.k5k2.kQk3.ki + 2k1.kQk2.k3ki.k5 
+ 2k1.kQk2.kik5.k3 — {d + 2)ki.kQk2.k5k3.ki] 

-\- 14 similar terms (4.3) 

where ki £ {k, I, q}. The full decomposition is clearly quite large. However, it is the appearance 
of the l/(d — 2) factor which is novel. In [17^118) the full set of three loop graphs in both massless 



and massive cases where a divergent F 



liver 
(3) 



(3) iiU' 



structure emerged, was noted. The sets of 



graphs appear to be the same. Though in [17] it is not fully clear which the actual ladder graphs 
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referred to are. However, the seemingly finite graph of Figure 4 was regarded as fully finite in 
all 7-string channels, |18]. In our present reconsideration it transpires that within the integral 
of the graph of Figure 4 there is a divergent contribution to the ^'(^^^ ^ r(3) ^i^^ channel but not 
for the r(o) <8' F^q) one. This derives from the pole — 2) in (j4.3p producing the massive Benz 
integral corresponding to the final graph of Figure 3. The key part is then 



(d — 2) Jkiq [k"^ — ra^] [/^ — m?] [q^ — m?] [{h — l)"^ — m?] [{k — q)"^ — m?] [(/ — q)"^ — m?] 

where the actual integral itself is finite in two dimensions. It remains after repeated application 
of (|3.22p and (j3.23p in the scalar integrals of (j4.3p . However, to have the complete divergence 
structure the integral needs to be evaluated since it will contribute to Z33. The remaining 
integrals with this — 2) pole in the F^^^^'^ '^{■i) ^vcr channel correspond to (j4.4p but with 
one or more propagators removed after application of ()3.22p and (I3.23p . These can be evaluated 
from the three loop techniques discussed earlier. In [I8j this contribution, (|4.4p . was overlooked 
since it was assumed that the parent integral with the internal momenta contracted was finite 
without noting the possibility of the l/{d — 2) factor deriving from the tensor decomposition. 
In relation to [ITJ we can only comment that in the massless version of (14. 4p the integral will be 
zero. However, given the totally different method of calculating the 4-point function of [18] in 
the massless contribution analogous to (14. 4p could possibly arise elsewhere. 




Figure 4: Three loop contribution to 4-point function. 



There remains the task now of evaluating the integral of (14. 4p . Although finite it clearly 
cannot be reduced to any of the three loop master vacuum bubbles already discussed even 
using, say, integration by parts. Instead we have had to resort to the more extensive experience 
of four dimensional vacuum bubble diagrams and specfically the Benz graphs discussed in [25]. 
To promote (j4.4l) to four dimensions we exploit Tarasov's observation of relating d-dimensional 
integrals to {d + 2)-dimensional ones, \2,Q\ I27j . Moreover, this is straightforward to do via the 
Tarcer package, [28], written in Mathematica for the basic two loop self energy topology 
given in Figure 5. Specifically one feature of Tarcer is that one can relate this two loop self 
energy graph in d-dimensions to that in [d + 2)-dimensions. This is a subgraph of Figure 4 with 
nullified external momenta and given that this is a three loop vacuum bubble, the final three 
loop integration measure can be rewritten as 



d'^k 



2iTd 



d'^+^k 1 
(27r)'^+2 



(4.5) 



in our conventions since the two loop subgraph will clearly be a function of k'^ only. From (j4.5p a 
massless propagator will appear in the higher dimensional integral. Since all the lines of Figures 
4 and 5 are massive and both final integrations involve functions of the square of the momentum, 
then we find the relation between the d-dimensional massive Benz graph and similar topologies 
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in two dimensions higher is 



Be(l, 1,1, 1,1,1, , , , , , , d) 



12m' 



-I I J\k^) 



+ 



TTd{d-l){d-2) 



'imp' Jk [fe^ — m'^] 

r\ rt rt rt r\ r\ 

Be(l, 1, 1, 1, 1, 1, m ,m ,m ,m ,m ,m ,d + 2) 
- Be(l, 1, 1, 1, 1, 1, 0, , m^,m^, m^,m^,d + 2) 



(4.6) 



where we define 



Be(a, /3, 7, p, A, m^, ma, mg, m^, m^,mQ, d) 

— .-3 



kiq - mf]"[/2 - mlf[q^ - mlY<[{k - If - ml]P[{k - qf - ml]\{l - qf - mjf 

(4.7) 

and emphasise that indicates a d-dimensional integration. The key part is the piece which 
represents the difference in two Benz topologies in (d + 2)-dimensions where one is completely 
massive and the other has one massless line. However, since these are multiplied by (d — 2) 
then in our e expansion relative to two dimensions we note that the leading term of each is 0(e) 
meaning that 



kiq [k"^ — m'^][P — m'^][q'^ — m'^][{k — If — m'^][{k — qf — m'^][{l — qf — m^] 
C(3) 



1927r3m6 



+ 0{e) 



(4.8) 



from (I3.18P and (IS.lOh . This is because whilst each of the two {d-\- 2)-dimensional integrals are 
divergent in four dimensions due to the presence of a simple pole in the regularizing parameter, 
the difference in (14. 6p is finite and the residue is independent of the masses in either Benz 
topology, [25] . 




Figure 5: Two loop self-energy topology central to Tarcer. 

With this observation all the ingredients are assembled to repeat the full three loop renor- 
malization of the 4-point function of (|2.ip . In [7J only the r(o) T(o) part was isolated but this 
was sufficient to deduce the /3-function at three loops. It is satisfying to record that we have 
verified the previous three loop MS result of [3 IS]. However, by contrast we find that a different 
renormalization constant from |17) and [18] emerges for Z33. We find 



^33 



C(3) _ J_ 
64 48 



whence 



"303) 
64 



1 " 

16 



(4.9) 



(4.10) 
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Though there is universal agreement on the rational part of (I4.10p , [17^ [T8] , only the contribution 
from the diagram of Figure 5 to the T'^^-j^ (^T (^^-^ channel produces the irrational piece thereby 
confirming the overall structure observed in [17]. However, rather than finding that we produce 
one of the previous values for Z33 we are in the seemingly unfortunate position of finding a new 
alternative. To determine which is correct and consistent it will transpire that the four loop 
mass anomalous dimension is the correct testing ground for this in the context of ()2.8p . 

5 Four loop vacuum bubbles. 

In this section we return to out initial aim and summarize the evaluation of the underlying four 
loop vacuum bubbles required for the mass anomalous dimension. For the 2-point function there 
are 18 distinct topologies and 36 Feynman diagrams to be considered. Of these topologies 14 
involve snail insertions in one way or another and hence their determination is in effect relegated 
to the straightforward extension of the three loop topology discussion. One effect of a snail is 
to produce two propagators on a line of a three loop graph but this can be reproduced by 
differentiating that line with respect to m^. This is also a reason why the three loop vacuum 
bubbles were required to be evaluated to the finite part exactly or left in terms of A(0) and 
other known integrals whose e expansion could be substituted when required, if at all. Several 
topologies contributing to the 2-point function, however, have a more demanding evaluation. 
These are illustrated in Figure 6 and we concentrate on these for the main part. Essentially the 
main complication now derives from rewriting the scalar products of internal momenta in terms 
of the propagator structure. For all the integrals which result we used several interconnected 
techniques. 




Figure 6: Several four loop diagrams contributing to the 2-point function. 

First was the use of the Tarcer package, ^28j, again, particularly for the third and fourth 
graphs of Figure 6. Clearly the last graph contains the two loop self-energy topology of Figure 5 
as a subgraph and the third has a similar two loop subgraph but with one line removed. Unlike 
the properties of Tarcer we described previously, the feature exploited in this instance is the 
ability to relate diagrams with different powers of the propagators in Figure 5 to that with unit 
power. Further, TARCER reduces integrals involving powers of the scalar products kl, kp and 
Ip where k and / are internal and p is the external momentum in Figure 5. The point is that 
the Lorentz tensor reduction for these situations can only be performed by this route. Any one 
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loop subgraph of Figure 5 will involve three external legs and the invariant decomposition in 
this case is too intricate. By contrast where possible we did exploit the Lorentz structure of 
subgraphs with one internal momentum flowing through it which can be regarded as a 2-point 
function external momentum for that subgraph. Then integrals can be rewritten using results 
such as 



1 



{d-l)^Ji [P - m'^][{k - If 



fc2 



fc2" 



.(5.1) 



The outcome of the Tarcer implementation is to reduce these more complicated tensor integrals 
to a set of master scalar four loop vacuum bubbles since the resulting combination of internal 
momenta allows for the repeated application of ()3.22p and (I3.23p . 

The use of (15. ip and Tarcer though may appear to introduce potential infrared difficulties. 
However, it transpires that in the full sum of all contributing pieces to a Feynman graph it can 
be checked that no integral retains an unprotected factor of l/k"^ which would give an infrared 
divergence upon integrating over the internal momentum k. For one instance checking this 
proved to be a tedious non-trivial exercise which we document for completeness. In all bar the 
second graph of Figure 6 the following combination of integrals emerge 



J{k^) - J(0) 



A(fc) 

A;2 



(5.2) 



Clearly each could be infrared divergent but the above combination always appears. Defining 

kii 



k [A;2 — m?Y[{k — pf — m?] 



then one can show 



1 



((i-4)J(p2 



1 r 



2' / / 2 
Using this and integration parts in (|5.3p one finds 



p^K{p^) + J(p2) _ j(o) 



Va = {d-3)i J^^J{k^)A{k^) 



2i / A{k)K{k' 
Jk 



ki P[{k-lf-m 



212 



(5.3) 



(5.4) 



(5.5) 



The final integral can be reduced using Tarcer if one regards the k momentum as external 
to the self-energy graph of Figure 5 and J{P) is replaced by the Feynman integral of (j3.5p . 
Consequently, Tarcer produces 



(d-2)2(d-4)l2 



{d-2){d-3) 



I /2[(A;-/)2 _m2]2 



1 



2(d-3)m2[A:2-m2]2 2m2[fc2-m2] Ji p[{k - I) 
^ (6^-4) (3d -8) .2 r 1 



212 



+ 



[k^ — m 
[{d-2)[k'^ - m 



iq [P — rrP][{k — (7)2 — rrp] 
21 8(d - 4)m2] 



[A;2 



m 



212 



X I 



1 



iq [/2 — rrPY[{k — qf — m? 



(5.6) 



The benefit of rearranging the two loop integral of the left hand side is to isolate the potential 
infrared singularity into a simple term on the right hand side. Moreover, the appearance of 
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powers of — m^] will lead to simplifications when substituted back into the expression 

for Va and the term with the singularity will actually combine with the first term of (j5.5|) to 
produce Va but with a factor of (d — 3). Hence, evaluating the remaining integrals of ()5.6p in 
the context of ()5.5p one arrives at the expression 



Fa = -if A{k)K{k'^ 
Jk 



(d-2)2/^A(0) 
2((i- 3)m2 



which has no potential infrared singular term. Moreover, each term of the right side of this 
is ultraviolet finite in two dimensions. So, in fact, when the combination Va appears in our 
computation, it can actually be dropped as there is no contribution to the renormalization of 
the mass at four loops. 

Having completed the tensor reduction of the scalar propagators, all that remains is the 
evaluation of a set of divergent four loop master integrals akin to those discussed earlier. Most 
of these are elementary given the results of (13. 7p . ()3.1ip and (j3.12p and the observation that 
as all integrals are infrared finite then one can ignore those four loop ones which are clearly 
ultraviolet finite by the usual counting rules. Though one integral is worth recording and that 
is 



' ' ' {2d -5) Jk ^ ' {2d-5){d-2) Jk 

because in the determination of this relation a pole in (d — 2) emerges in the standard d- 
dimensional manipulations such as differentiating the original integral with respect to and 
using (j3.7p . This pole gives rise to a problem similar to that discussed for (j4.4p but with a 
simpler resolution since we merely apply the technique used to deduce (|3.18p and (j3.19p , to find 

The need to evaluate i Jf^ J'^{k'^) to the finite part as well is also illustrated by this equation. 

This completes the discussion of the construction of the relevant basic Feynman integrals. 
For each of the eighteen topologies a Form module was created within which the algorithm 
to break the original Feynman graphs up into its basic components was encoded. The tedious 
identification with the above results together with the remaining more elementary ones were 
also contained in each module. Finally, prior to summing all the results from the 36 four loop 
diagrams, the e expansion of / and other integrals were evaluated to the appropriate order in e. 
The resulting sum produced the divergent part of the mass component of the 2-point function 
to the simple pole in e as a function of the bare parameters. 



6 Four loop renormalization. 

The final piece of the calculation rests in determining the overall renormalization constant Z^a at 
four loops in MS. However, prior to this we must consider the full theory. To this point we have 
tacitly assumed that only the basic 4-point vertex of (j2.ip is responsible for all the Feynman 
diagrams we have discussed. The presence of the generated evanescent operator in ()2.9p needs 
to be included. As noted earlier since the operator appears with a coupling the effect of this 
operator cannot arise before four loops. Therefore, we now have to include the additional graph 
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of Figure 7 where the chxle with a cross in it denotes the insertion of the operator O3 with its 
associated renormahzation constant Z33. The integration routine to determine its contribution 
is the same as that for the original vertex except that one has to first replace the r^^^*^ matrices 
by the corresponding string of ordinary 7-matrices. 



Figure 7: Contribution from evanescent operator to the four loop mass renormahzation. 

With this additional graph included the overall renormahzation constant is extracted using 
the standard method for automatic Feynman diagram computations developed in [29]. Briefly 
one computes the Green's function of interest as a function of all the bare parameters such 
as the coupling constant and the mass. Then the renormalized parameters are introduced by 
the rescaling defined by the renormahzation constants. In the present context these are the 
renormahzation constants leading to the naive anomalous dimensions as defined in (j2.9p and 
()2.10p . This rescaling in effect reproduces the counterterms to remove subgraph divergences. 
Moreover the Green's function is multiplied by the associated renormahzation constant which 
in our case and conventions is Z^Z^- As the former is already known, [9], then the divergences 
which remain in the 2-point function are absorbed by the unknown pieces of Zm- We recall that 
at four loops the anomalous dimension of ^ corresponds to the naive anomalous dimension 
■y{g) since there is no contribution from the graph of Figure 7 in the wave function channel. 
Therefore, having followed this procedure we find the naive mass anomalous dimension in MS is 



At this stage several comments are necessary. First, there are several checks on the underlying 
renormahzation constant itself. Whilst the evanescent operator issue arises at four loops, it will 
manifest itself in the simple pole in e of Zm- Therefore, the quartic, triple and double poles 
in e are in fact already predetermined by the structure of previous loop order poles from the 
renormahzation group equation. For (|6.ip we have verified that this is in fact correct. One other 
useful check was the explicit cancellation of divergences of the form A(0)/e" for n = 1 and 2 at 
four loops. This is non-trivial since, for instance, A(0) arises at three loops both associated with 
a simple pole in e and in the finite part. Therefore, one needs to write A(0) as a formal expansion 
in powers of e prior to the rescaling of the bare quantities. This is because the 0(1) piece at three 
loops will be multiplied by 1/e poles. Moreover, since A(0) has dependence (m^)'^~^, then this 
has to be explicitly factored off since this mass is bare and needs to be renormalized too. Once 
written in this way we have checked that the poles in e involving the 0(1) and 0(e) residues 
stemming from the e expansion of A(0) do indeed cancel completely. 

Again one can partially check part of (j6.ip from another point of view. In [9] the structure 
of the mass anomalous dimension has been given in the large expansion to 0(1/A^^) based on 
the results from a series of articles [301 EH 1321 ESI EH ES ES] . Again at this level of expansion the 



Q 



7m (5) 



-(2iV-l)^ + {2N-l)^ + (4iV-3)(2iV-l)-|-3 

+ \{48N^ - 384A^2 + 492A^ - 138)C(3) - AON^ - 72N'^ + WON - 81 
+ 0(/) . 
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evanescent operator is not manifested and so the 0{1/N'^) coefficients of the mass anomalous 
dimension which are given there at four and higher loops in fact equate to those of the naive mass 
anomalous dimension ^m{g)- In other words if it were possible to compute the critical exponent 
corresponding to the mass anomalous dimension at the d-dimensional fixed point of the theory 
at the next order in large A^, 0(1/A^^), then unless the effect of the evanescent operator could 
be included, it would not correspond to the true mass anomalous dimension, [9]. From the 
expression given in [9] we note that when the same convention is used, that part of (16. ip at four 
loops which corresponds to the 0{1/N'^) piece agrees precisely with [9]. This is a reassuring 
cross-check on a significant part of our four loop computation since, within the computer setup, 
one can examine the A^-dependence multiplying all the basic integrals which we have had to 
compute for all topologies. The vast majority are at least touched by a quadratic or cubic in 
which are related respectively to the 0{1/N'^) or 0{1/N) large N piece already determined 
in [9]. For the small number of remaining pieces which have linear factors in N we have been 
careful in evaluating the corresponding, though invariably simple, vacuum bubbles. Therefore, 
we are confident that (16.11) is correct. 



One clear problem remains which is related to the structure of the expression (j6.ip . Unlike 
the previous orders the four loop part does not vanish when N = \ which corresponds to the 
free field theory. Moreover, it transpires that of the eighteen underlying topologies only the 
graphs for one do not vanish for this value for N . (Though actually the parts from the second 
and third graphs of Figure 6 cancel between each other which is similar to what occurs at three 
loops for analogous graphs.) The topology which gives a contribution for = ^ is the final 
graph of Figure 6. However, given our discussion in several places concerning the evanescent 
operator, the resolution is clearly straightforward. More concretely one can see the evidence for 
this if one evaluates ()6.ip at A^ = ^ to find 



imig) 



This is the piece which needs to be cancelled in order to have a mass dimension consistent with 
a free field theory. Indeed this is the relative combination of rationals and C(3) which our three 
loop 4-point function renormalization reevaluation produced. Therefore, using (12. 8p and (I4.10p 
we can derive the true mass anomalous dimension as 

2 3 

^m{9) = -{2N-l)^ + {2N-l)^ + (4Ar-3)(2Ar-l^ ^ 



27r ^ ^87r2 ^ ^327r3 

2 .... , .... a' 



12{2N - 13)(A^ - 1)C(3) - - A<dN + 57 (2A^ - 1) 



+ .... , „,j 334^4 

+ 0(5') • (6.3) 

Clearly this has the correct expected A^ = | property and given our earlier checks on (|6.ip we will 
regard (16. 3p as the completion of our original aim. Also, it is worth stressing that the discrepancy 
in the 4-point function renormalization has now been crucially resolved simultaneously. It turns 
out that neither of the previous expressions for (33(g), [iTl |18], could be correct to preserve 
the vanishing of ^m{g) in tbe free field case. So we can regard this mass anomalous dimension 
calculation as also a non-trivial check on the full three loop MS renormalization. 



7 Discussion. 

We have completed the four loop renormalization of the mass anomalous dimension of the Gross- 
Neveu model in the MS scheme. Despite the lack of multiplicative renormalizability when the 
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Lagrangian is regularized dimensionally, it has been possible to compute an expression which 
passes all possible internal checks. Not least of these is the correct implementation of the 
projection formula formalism of [131 [H] which has been justified by the consistency with the 
free field case. Concerning this the previous attempts to deduce Psig) appear to indicate that 
the only approach which is truly reliable for renormalizing the model is the one where there is a 
non-zero mass. This seems to be the conclusion one must draw from the origin of the necessary 
C(3) part missing from [17] required to balance the discrepancy of ()6.2I) . Given these remarks 
one possible extension would now be to repeat the derivation of the mass anomalous dimension 
at four loops in other 4-fermi models in two dimensions. Whilst considering the most general 
possible interactions involving Oq, Oi and O2 of [131111] would perhaps be too ambitious, there 
is the interesting case of the non-abelian Thirring model, [371 EH]- The seed interaction involves 
Oi but includes colour group generators too. It has been renormalized at three loops in MS in 
|18| and the four loop wave function is also known, |10j . Though in light of our comments on the 
4-point function in the Gross-Neveu model, the corresponding 4-point function renormalization 
would clearly need to be reconsidered to deduce the correct evanescent operator /^-functions. 
One motivation for determining the mass anomalous dimension in the non-abelian Thirring 
model would be to examine the colour group Casimir structure of the final expression since, 
given the similarity with QCD, it is thought that it should involve the same structures as the 
corresponding expression for the quark mass anomalous dimension, |391 140j. This was the case 
for the wave function, jlOj . 
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